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1. INTRODUCTION. Mathematical studies of the dynamic stability of 
aircraft began essentially with Bryan's "Stability in Aviation" pub- 
lished in 1911 [1]. This analysis appeared at the very beginning of 
heavier- than-air flight itself, and has remained the foundation for 
practically all subsequent studies of the subject. B. Melville Jones 
[2] reported on progress 25 years later and succinctly stated the sub- 
ject's principal task: 

"Given the shape of the aeroplane and the properties of the 
air through which it moves, the air reactions X, Y, Z, L, 

M, N depend on the motion of the aeroplane relative to the 
air; that is to say upon the six variables U, V, W, P, Q, R 
and their rates of change with respect to time. In prac- 
tice, the principal difficulty lies in determining the 
relationships between X, Y, . . . and U, V, . . ." 

The establishment of these relationships with sufficient realism is 
what we now recognize as the province of mathematical modeling. 

Bryan's formulation, which originated the subject, has at its core 
the assumption that the aerodynamic forces and moments developed at a 
given instant are functions only of the instantaneous values of the 
variables that determine forces and moments in a steady flow. When in 
addition a linear dependence of the forces and moments on these 
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variables is assumed, the equations governing the motions of aircraft 
reduce to a set of ordinary differential equations having constant 
coefficients. With the form of the equations thus established, the 
mathematical modeling problem is completed, and the study of stability 
becomes synonymous with the study of the coefficients. The nature and 
the determination of these coefficients, the "stability derivatives," 
have been the central concerns of experimenters and analysts alike 
through the ensuing years. 

Later investigations into the transient behavior of the forces and 
moments in response to sudden changes in the flow around the aircraft 
led researchers to recognize that the forces and moments at an instant 
were dependent not only on the instantaneous values of the flow varia- 
bles but also on their past values (cf. [3] for a comprehensive summary 
and bibliography). The concept of transient aerodynamic force and 
moment responses to step changes in the flow variables, i.e., of 
"indicial functions," coupled with the notion of superposition, led to 
a new formulation of the equations of motion [4]. This formulation is 
exact in principle within the assumption of linearity and the equations 
of motion take the form of integro-dif ferential equations. However, 
reduction of the equations to equations correct to within a first-order 
dependence on time-rates-of-change of the variables restores the form 
of the original ordinary differential equations, which now include terms 
that account for the past within the order of the approximation [5). 

From the standpoint of mathematical modeling, exploiting the concept 
of a linear indicial response was an important step, since it described 
the process of generating an aerodynamic response to an arbitrary 


motion from a few aerodynamic indicial responses by using superposition 
integrals. It is from this idea that the modeling approach derives its 
economy in the treatment of time-history effects. Further, use of the 
superposition integral overcame in a concise way the first and most 
objectionable assumption of Bryan’s formulation: the aerodynamic forces 

and moments could now depend on the history of the motion. Mathemati- 
cally speaking, expressing the aerodynamic forces and moments in the 
form of linear superposition integrals replaced Bryan’s functions by 
linear functionals . 

As time passed, the continual expansion of the aircraft performance 
envelope brought new aerodynamic phenomena, such as shock waves and 
concentrated vortex flows, into play. It became apparent that the ways 
in which these phenomena influenced aircraft motion could not all be 
encompassed by even the exact linear formulation of the aerodynamic 
forces and moments. A reformulation of the aerodynamic system that was 
free of the remaining assumption basic to the original formulation was 
now in order. It could be anticipated that the reformulation would 
result in the replacement of the superposition integrals, which were 
linear functionals, by suitable nonlinear functionals. As it turned 
out, this task could be accomplished easily by the adoption of Volterra’s 
original conception of a functional [6]. Functional analysis was used 
to construct a framework within which the indicial function could be 
reformulated as a nonlinear functional [7]; the result was a new defini- 
tion for the indicial function that did not depend on a linearity 
assumption. This definition led naturally to the derivation of integral 
forms for the aerodynamic forces and moments which were the anticip ted 



generalizations of the superposition integrals. As in the case for the 
exact linear formulation, it was found possible to reduce the results 
to simpler, more practicable forms by virtue of the low angular rates- 
of-change characteristic of aircraft motions. Subsequent papers [8-16] 
revealed the implJ cations of these results in regard to general motions, 
experiments, and a variety of modeling questions. A comprehensive 
summary of the work to date is available [17]. 

In recent years, we have become increasingly aware that our extensive 
use of concepts from functional analysis in the modeling of systems 
governing the dynamical motions of aircraft has been matched or sur- 
passed by workers in a wide variety of fields. Noteworthy in this par- 
ticular setting has been the body of work which has developed in several 
branches of electrical engineering (for a comprehensive survey, see the 
entire contents of [18,19]). We consider it important that the common 
features of the efforts in these various fields be brought to the fore 
for everyone’s benefit. To that end, we shall take this opportunity to 
recast the basic ideas of the approach we have taken in the modeling of 
nonlinear aerodynamic responses in a way that we hope will be compatible 
with some of the approaches taken in electrical engineering. We shall 
emphasize the physical aspects of our approach to clarify its relation, 
in particular, to the body of ideas underlying the use of nonlinear 
functional expansions (cf. [20] for an excellent exposition of the latter 
work). Finally, we shall try to show how our analysis can be extended 
through its natural connection with ideas from bifurcation theory. 



2. MODELING INFLUENCE OF PAST MOTION BY PULSES. For clarity of 
exposition, it is convenient to adopt the two-dimensional wing as illus- 
tration. Results will have more general bearing, however. 

Let the wing move away from a coordinate system whose origin is fixed 
in space at the center of gravity at a time £ - 0. The distance 
traveled by the center of gravity along the flightpath is measured by a 
coordinate s. Let the center of gravity move at constant velocity V 0 , 
so that the trace of its path, plotted against time £, is a straight 
line. This is shown on Fig. 1. The wing is allowed to undergo changes 
only in the angle of attack a, where a is the angle between the 
velocity vector and the wing chord line. Projections of the leading 
and trailing edges of the wing onto the plane containing the velocity 
vector are maximum when a - 0. These maximum projections also trace 
out straight lines on Fig. 1, parallel to the trace of the center of 
gravity. 

As illustrated in Fig. 1, let the angle of attack a be zero for all 
time £ except at £ = £ 1# where a pulse occurs of amplitude a(£ a ) 
and infinitesimal width A£ x . Consider a measuring point s on the 
wing at a time t subsequent to £ x . The loading at the point is 
influenced by all disturbances that originate in the past and are able 
to reach the point at the same time t. Each disturbance is propagated 
at the local speed of sound, and hence, in a plot such as Fig. 1, the 
zone of its influence is bounded by projections of the rays of an 
approximately conic surface whose origin is the point of the disturbance. 
Only disturbances whose cones include the point (s,t) in question can 
influence the loading at the point. Thus, a certain conic surface, 
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directed backward in time from the point (s,t) will include within it 
all points in past time whose disturbances are able to influence the 
loading at s at time t. The projection of such a conic surface is 
shown in Fig. 1, where, in the present case, the only disturbances that 
exist are those originating from the pulse at £ - It will be seen 

that only disturbances originating from the shaded part of the pulse 
can influence the measuring point. Also note that if the elapsed time 
t - between pulse and measuring point is held fixed and if the same 
pulse and the measuring point are translated together to new positions 
on the wing in such a way that the trace of the measuring point remains 
parallel to that of the center of gravity, then nothing changes in the 
form of the loading at the measuring point. This behavior is captured 
by writing the loading at the measuring point in the form of a Taylor 
series: 

AC p (s,t) = AC p (x,t - C 1 .a(C 1 )) dir * ^ a n ( X,t ' ?i)t3(f; i )] n AC 1 

n ( 2.1 

where loading AC p is the difference between pressure coefficients on 
the lower and upper surfaces, and x is the distance of the measuring 
point from the maximum projection of the leading edge (Fig. 1). It will 
be noted that, as required, AC p remains constant when x, t - C l , and 
a(C x ) are held fixed. The second of the forms in (2.1) will be used in 
the subsequent analysis to distinguish between direct (subscript dir) 
and interference (subscript int) effects. Equation (2.1) holds under 
the principal assumption that at least a limited range of a will 
exist in which the loading will depend analytically on o. Further, in 
more general circumstances, such as in accelerated motion where the 
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trace of the center of gravity will be a curved line, the dependence on 

elapsed time t - alone will not hold; the a will then depend on 

1 n 

t and separately. 

Now let us consider the response at the measuring point to a pair of 
pulses located at and £ 2 with < t. Here, in addition to the 

direct influence of each of the pulses acting as if in isolation, the 
interference between the pulses will also influence the loading. The 
interference effect can be written in a form resembling a product of 
responses to single pulses 


AC 

P 


(s,t) 


int *2 



m,n 


- 5 x .t - C 2 )[a(C 1 )] m [a(C 2 )] n A? 1 A£ 2 

( 2 . 2 ) 


where the subscript (int ,2) is meant to be read as '’interference between 
a pair of pulses.” With the addition of the direct influence of the 
two pulses, the loading at (s,t) takes ths form 


AC (s,t) - AC (x,t - 5 1> a(C 1 )) dlr + AC (x,t - Z 2 >*<Z 2 )) dir 


* a y s,t) lnt,! 


(2.3) 


The process of adding pulses can be continued indefinitely in the 
same way. At the next stage, e.g., the Interference between triplets of 
pulses must be considered as well as that between pairs. Going to the 
limit of a continuous distribution of pulses starting at time f‘ * 0 
yields a summation of multiple integrals having the form 


AC p (s,t) 


AC + AC 

^dir ^int,2 


+ AC 

P 


int , s 


+ • • 


(2.4) 
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with 


AC 

^dir 



a n (x,t - C 2 ) [a(C 1 ) ] n d^ 


(2.5) 


AC 

P int ,2 



[a(5 2 )] m dC 2 




f n ,t - 5 2 )[o U 1 )] n dE, 1 

u.o) 


AC 

P int , 3 



in , n , p o 


[ct(C 3 )] m dC 3 


I [«(C 2 )] n dC 2 


X I C mnp^ X,t - f 'i’ C * *1* ~ e s >l«CC x >I P dC 2 (2.7) 

*0 

Written as a nonlinear functional expansion, (2.4) represents the load- 
ing at the point (s,t) in response to an arbitrary variation of a 
over the time interval zero to t. The form of (2.4) confirms an 
important point made in [20]. It will be noted that a summation of the 
leading terms in (2.5,2 .6 , . . . ) forms a Volterra series [6,20]. The fact 
that there are terms left over confirms that the k priori adoption of a 
Volterra series to represent the loading would have been insufficiently 
general to accommodate the Taylor series form of the dependence on angle 
of attack. 


3. FORMATION OF INDICIAL RESPONSE. Given (2.4), one can now use it 
to form the indicial response as we have defined it in [7]. To indicate 
the form of the result, it will be sufficient to consider terms in (2.4) 
only through the series representing AC p ^ t 2 * Two morions need to be 
considered. In the first, the wing undergoes the motion under study 
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a(4) from time zero up to a time 4 ■ t, where x < t. Subsequent to 
i, a is held constant at cx(t). Thus, in (2. 4, 2. 5, 2. 6) 
a(0 - <*(0 ; 0 < 4 < t 

■ a(r) ; 4 > t 

The direct and interference contributions take the forms 


AC 


'dir 


AC 


int i 


m /\ ® n < x,t " 
n *o 

+ ^ j a(t) ] n f a n ( x » c ~ £ x )d4j 

VT T 

2 j f a ^2>J md ^2 J b mn [o(C i ),ndC x 

* fflTn o o 

*WJ 

mVn t o 

f “■ J' 


b »n 


(3.1) 


(3.2) 


In the second motion, the wing undergoes the same angle of attack 
history a(4) up to time x. Subsequent to t, the angle of attack is 
again held constant, but is given an incremental step change Aa over 
its previous value of <x(t). Thus, in the second motion, 
a (O ■ a(C) ; 0 

- u(t) 4- Aa ; 4 


< C < t | 

> t J 


(3.3) 



i 
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The direct and Interference contributions become 


a n [a ( c 1 >] n df ; 1 + ^° 


(t) + Aa] 


a n d ^ 




(? 2 )] m dC 2 


b [a^)] 

mn 1 i 


m , n o 


•Vmt) + Aa] m J d «2 J b mn [a( 5 i )]n dC i 

m,n t o 


+ 7 [«(t) + Aa] 


r-r 


b dr 
mn i 


The difference between the two loadings, with terms retained only of 
O(Aa), yields the incremental change in loading in response to the 
incremental step change in angle of attack: 


AC - AC 
P 2 


— ■ N l 


(x) ] n " A a n d^ 


^m[a(,)r‘J It, J 

m,n t c 

^ (" + nMaCT )]"*"' 1 J dC 2 J b fln dl 


b mn df (3.5) 

fWCX 1 


Equation (3.5) reveals the form of the indiclal loading response to a 
step change in angle of attack in terms of functional expansions. It 
will be seen that the first and third terms on the right-hand side of 
(3.5) do not depend on the past motion but only on the level of the 
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angle of attack a(x) at which the step was made. The second term 

depends on the past motion, however, since a(£ x ) with 0 < < t 

appears within the integral. The leading term of this past dependence 

has the form J d£ 2 J a(C 1 )b 11 (x,t - S x ,t - £ )d£ t . Dependence on 

the past thus arises from interference effects between pulses prior .o 

i, the origin of the step, and perturbation pulses of 0(An' originating 

subsequent to x. In the general nonlinear case, then, and just as 

before [7], the indicial response is itself a functional. In the limit 

as Act -+ 0, then (3.3) can be cast in the functional form 

AC - AC 
P P 

lim (s,t) - AC [a(C);x,t,x] (3.6) 

Aar*o P P a 

where it will be noted that, within the functional, dependence on t 

and x is indicated separately rather than on elapsed time t - x alone. 

It can be easily verified that the first and third terms in (3.5), which 

depend only on the level ol(x), indicate a dependence on t - x alone; 

however, as a consequence of its dependence on the past motion, the 

second term cannot be cast as a function of t - x alone. 


GENERALIZED SUPERPOSITION INTEGRAL. Just as before, (3.6) can 
be used to form a generalized superposition Integral for the response in 
ACp to an arbitrary angle-of-attack variation. The result is 


AC (s , t ) 
P 


- Aw p (s,0) 



‘iC [a(0;x,t,T] 
p a 


(4.1) 


By substituting (3.5) for AC [a(0 »x*t ,t J in the integral term in 

"a 

(4.1) and carrying out the integration, one will /erify that the form of 
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(4.1) is restored through terms of the series representing AC 

P int ,2 

This should suffice to demonstrate that our use of indicial responses 
and generalized superposition integrals is in fact compatible with the 
approach based on nonlinear functional expansions. 

5. PRINCIPAL SIMPLIFICATION. The simplification most instrumental 
in making our use of nonlinear functionals in the analysis of aircraft 
dynamics practicable has proved to be the reduction of the general 
aerodynamic force and moment responses to forms correct to within a 
first-order dependence on angular rates. This approximation is justi- 
fied in its application to studies of aircraft dynamics in view of the 
generally low reduced frequencies characteristic of aircraft motions. 

In the following, the steps involved in the reduction are reviewed in 
order to highlight the several advantages of the reduction as well as to 
set the stage for a discussion of the amendments required to accommodate 
the occurrence of bifurcation phenomena. 

Integrating the expression for loading in (4.1) across the chord 
eliminates the dependence on x and yields a form for the response In 
lift coefficient to an arbitrary angle-of-attack variation 

C L (t) - C L (0) + I C L [t,T;a(0] dx (5.1) 

i a 

Note that the order of the dependencies has been reversed within the 

functional to anticipate the enhanced role played by t and t in 

a 

the ensuing analysis. 

Consider first the behavior of the functional as the elapsed time 
t - x increases. It is clear physically that, with Increasing t - t, 
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the dependence of on the past motion a(£); ^ < t must fade away, 

a 

Thus, as t - t + 00 , approaches a function which is dependent only 

a 

on t * t and a(i), the level of the angle of attack at which the step 

was made. Additionally, we assume here that as t - x 00 , will 

a 

approach a unique, constant value corresponding to the lift-curve slope 
which would be measured in a steady flow. More precisely, we assume 
that fluctuations, which always exist, are small enough to be neglected 
in comparison to mean values resulting from ensemble averaging. It is 
this assumption that will require amendment in the consideration of 
bifurcation phenomena. Here, however, we invoke it explicitly by making 
the substitution 

C L ft ,t ;ci(0 ] * C L (°°;a(i )) - F[t,i;a(£)] (5.2) 

a a 

The quantity F is called the deficiency function (al) . Provided the 

steady-state term (°°;ct(T)) is in fact constant, F must approach 

a 

zero as t - t + °°. When (5.2) is substituted in (5.1), the product of 

(®;a(t)) and a(x) forms a perfect differential which can be inte- 

a 

grated, yielding 

C L (t) = C L (®;u(t)) - J F[t ,t ; ot (C ) ] dt (5.3) 

o 

where C^O;a(t)) represents the (mean) lift coefficient that would be 
measured in a steady flow with a held fixed at the instantaneous value 
u(t). Now we examine the dependence of F on the past motion a (5)* 

If a(£) can be considered an analytic function in a neighborhood of 
£ =* t (corresponding to the most recent past for an indicial response 
with origin at £ * x), then its history can be reconstructed, in 
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principle, from a knowledge of all of the coefficients of its Taylor 
series expansion about £ * t. Since a (5) is equally represented by 
the coefficients of its expansion, it follows that the functional F, 
with its dependence on a(£)> can be replaced without approximation by 
a function which depends on all of the coefficients of the expansion of 
a (5) at £ * t. Thus, F can be expressed as 

F[t , t ; a(£) ] = F(t , t ; a(i) ,a(x) , . . .) (5.4) 

Assuming now that a(x) is potentially large but that the rates 
d(x),d(x), . . . are all always small, we are permitted to expand (5.4) 
around the zero values of the rates, so that 

F [t ,t ; a(£) ] = F(t ,t ; a(x) ,0,0, . . .) + a(x)F^(t ,x ,a(x) ,0,0, . . .) 

+ «_Ill F^(t,T;a(T) ,0,0, 

+ a(x)F-(t ,x ;a(x) ,0,0, (5.5) 

Returning to (5.3) and retaining terras only to within a linear dependence 
on a(x), we get for the integral term 

I(t) * J F(t ,x ;a(x) ,0,0, . . .) dx (5.6) 

"o 

To the order of the approximation, only the first term of the expansion 
of F survives in (5.6). As a consequence, and as the notation shows, 
as far as the functional F is concerned, the past motion is simply the 
constant motion o(C) » a(x). On this basis, it is consistent now to 
say as well that F will depend only on elapsed time t - x, rather 
than on t and t separately. Expansion of a(i) about x ■ t and the 


dependence of F 
dence on a(t). 


on 


a(x) about a(t) yield, to within a linear depen- 
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I(t) * d(t) 


<x 

I 


F(u;a(t))du 


(5.7) 


where we have assumed F(u;a(t)) a 0 for all u = t - x > t a and t 

large enough so that t > t . Substituting this result for the integral 

€% 

term in (5.3) yields as the formulation for the lift coefficient 


C L (t) = C L ( ~ ;0t(t)) + i(t) V C L. (a(t)) 

o a 


(5.8) 


where 


C L (a(t))=-^f F(u;a(t))du (5.9) 

* J o 

and 1 Is a characteristic length (e.g. , chord length). 

The simple form of (5.8) contains a number of important advantages. 
Since C depends only on a(t), within the approximation the response 

L • 

a 

to any of a class of sufficiently slowly varying motions which arrive at 
the same value of a at time t will yield the same value of C . 

Ij • 

a 

In particular, a harmonic motion of infinitesimal amplitude about a mean 

value of a equal to the instantaneous value a(t) suffices to obtain 

C . This motion is the preferable one from both the experimental and 
a 

the computational standpoint. On the other hand, a finite amplitude 

harmonic oscillation about a mean value of a, say a , can be used as 

m 

well to obtain C. at the particular points in the motion where the 
L* 
a 

instantaneous amplitude a(t) is such that a(t) ■ a(t) - a^. This fact 
enables one to deduce the results for finite-amplitude oscillatory 
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motions from those of infinitesimal-amplitude oscillations, a consider- 
able advantage from the computational standpoint. Finally, and most 
important, as will be seen in the contribution to this collection by 
Hui and Tobak, the simplification provides the rationale that allows 
decoupling flow-field equations from the inertial equations of motion 
and using results for the aerodynamic terms in the mathematical model 
from known harmonic motions. 

6. CONNECTION WITH BIFURCATION THEORY. We intend to show that our 
use of indicial responses as a basis for arriving at the form of the 
aerodynamic force and moment leads naturally to the consideration of 
bifurcation phenomena. To fix ideas, we restrict attention initially to 
laminar flows and assume that the flows are governed by the Navier- 
Stokes equations. Further, we assume that the maneuvers and measure- 
ments can be carried out with sufficient precision so that they are 
effectively er\ or less and repeatable, eliminating the necessity of 
accounting for the presence of random fluctuations. The extent to which 
the analysis must be deepened to accommodate turbulent flows and random 
fluctuations will be addressed later. 

We begin by reconsidering the maneuvers required to form an indicial 
response. As before, two maneuvers are involved, both beginning at 
£ * 0, z.id constrained at £ - t, and differing only in the step imposed 
on t~.e second maneuver at £ * x. For each maneuver, the lift coefficient 
is measured at a time t subsequent to T. If we assume that the dif- 
ference between measurements at time t, AC^(t), divided by the magnitude 
of the tftep, An, exists and is unique in the limit as Act 0 for all 



values of £ ® t > t » then we define this limit as the indicial response 
in lift coefficient per unit step change in a. Hence, the basic 
assumption on which the definition of the indicial response rests will 
fail when the indicial response ceases to exist or to be unique. A very 
natural way of invalidating the assumption is through the mechanism of 
instability, which we have not considered in this context before. It 
is here that the possibility exists of extending our analysis by incor- 
porating ideas from bifurcation theory. 

Consider the first of the two maneuvers involved in the formation of 
the indicial response. The angle of attack attains a constant value 
ol(t) subsequent to t and it is reasonable to expect that the flow 
field at the subsequent time t will approach an equilibrium state that 
corresponds to this fixed boundary condition, as the elapsed time 
t - t •> ot> # in all of our previous analyses, we have assumed that as the 
flow field approached the equilibrium state it became time-invariant, 
which meant that the corresponding lift coefficient C^(t) approached a 
unique constant value C (°°;a(T)) as t - t ■+ 00 . As long as this was 

L 

true, it was reasonable likewise to expect that an incremental change 
in a(r) of 0(Aa) would result in an incremental change in C^(»;a(t)), 
likewise of 0(Aa). We now recognize that this will be true as long as 
the time-invariant equilibrium state represented by C^OsaCr)) is 
(asymptotically) stable to small perturbations in a. It can happen, 
however, that as a(x) is increased in small Increments, a critical 
value of o(t) can be reached at which the stationary equilibrium state 
represented by C^(<»;a(T)) will no longer be stable to small perturba- 


tions in a. Of the new equilibrium states that are possible, the 
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system will seek one which can remain stable to small perturbations. 

Now this is precisely the situation that bifurcation theory is designed 
to address. Bifurcation theory tries to classify and characterize the 
properties of the new equilibrium states that can arise when the given 
equilibrium state becomes unstable. Of the types of bifurcation phe- 
nomena that are possible, perhaps the most typical in aerodynamics is 
the "Hopf" type, which is characterized as follows: A previously stable 

time-invariant equilibrium state is replaced by a time-varying oscilla- 
tory equilibrium state. Physically, the usual origin of such a large- 
scale oscillatory state is the onset of vortex-shedding. Of the many 
examples, we cite here stall on airfoils when the angle of attack 
exceeds a critical angle [21,22] and the wake of the flow past a cylin- 
der when the Reynolds number exceeds 50 [23]. Bifurcation theory gives 
us the means to incorporate these phenomena within a rational frame- 
work, consequently, the possibility of accounting for those critical 
points in maneuvers where sudden and dramatic changes in flow structure 
may occur. 

To conclude, we indicate our current thoughts on the directions in 
which the analysis must be extended to acknowledge the important effects 
of random fluctuations and turbulent flows. The issue has separate 
experimental and computational components. 

From the experimental standpoint, the presence of random fluctuations 
in the maneuvers and the measurements is a practical question which 
arises even with strictly laminar flows. As we have already indicated, 
however, (cf . , in particular, [12]) adoption of ensemble-averaging allows 
us, in principle, to acknowledge the presence of random fluctuations 


without otherwise having to alter the analysis. On this basis, one can 
recognize the probable existence of bifurcation phenomena by the exami- 
nation of the mean value of the lift coefficient in the equilibrium 
state, say C^(°°;a(T)). Values of a at which C l («;u(t)) is double- 
valued, or at which the slope of C L 0»;a(x)) with a is discontinuous, 
are the signs of probable bifurcation phenomena. We have already noted 
the existence of the latter symptoms in several of our previous studies 
(cf. [15,16,17]) and, in [17], we have devised a special scheme for 
treating them in the particular instance where they reflect the presence 
of hysteresis in the equilibrium flow. We now recognize the possibility 
of incorporating such special treatments within the more general frame- 
work that an analysis based on ideas from bifurcation theory will provide. 

Consideration of the question of random fluctuations and turbulent 
flows from the computational standpoint brings to the fore an additional 
issue. While laminar flows could be said to be governed by the Navier- 
Stokes equations, even in the presence of random fluctuations, account- 
ing for the presence of turbulent flow in computations centers around 
the problem of having to define the equations governing the flow. This 
is the turbulence modeling problem. In current practice, it is gener- 
ally agreed that any particular realization of a turbulent flow could 
be modeled, in principle with sufficient accuracy, by a solution of the 
Navier-Stokes equations. However, the existence of rapid and apparently 
random fluctuations in the flow makes it mandatory that the equations 
be averaged to suppress the appearance of the rapid fluctuations. The 
averaged equations, called Reynolds-averaged Navier-Stokes equations, 
are then closed by the installation of a suitably chosen closure model. 
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This effectively casts the equations in a form similar to that of the 
original Navier-Stokes equations. Therefore, in application to modeling 
the equilibrium flows associated with the formation of indicial responses, 
we should expect the modeled turbulent flow equations to mirror behavior 
previously captured by the Navier-Stokes equations in application to 
laminar flows. Thus, typically, when they are applied at low values of 
a(x), modeled turbulent flow equations should yield solutions for equi- 
librium flows that are invariant with time. But, just as before, time- 
invariant equilibrium flow solutions that had been stable for values of 
ot( t) below a critical value should be expected to become unstable upon 
exceeding the critical value, and to seek a new branch of stable solu- 
tions. One possibility is a branch consisting of time-varying oscilla- 
tory solutions. Encouraging evidence is available that modeled turbulent 
flow equations in fact can be sufficiently general to exhibit such 
instability and bifurcation phenomena. In particular, the results of 
Levy ([24-26]; see also the discussion in [27]) for transonic flow past 
a biconvex airfoil show the typical Hopf-type bifurcation that reflects 
the onset of vortex-shedding in the wake, when either Reynolds number, 

Mach number, or lift coefficient exceed critical values. An analogous 
problem involving the occurrence of aileron buzz at transonic speeds 
when either Mach number or angle of attack exceeds critical values, has 
been treated by Steger and Bailey in [28]. 

The concerns of turbulence modeling research and those of research in 
the modeling of aerodynamic responses converge on the issue of the 
bifurcation behavior of equilibrium flows. In judging the importance of 
this issue, one should note that bifurcation phenomena reflect the 
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occurrence of sudden and potentially dangerous changes in flow structure 
during dynamic maneuvers. A resolution of the issue hinges in great 
part on the correct identification of these occurrences as the relevant 
parameters (e.g., angle of attack, sideslip angle, Reynolds number, Mach 
number) range over their respective envelopes, and their capture by 
means of computations based on modeled equations of turbulent flow. In 
application to two-dimensional flows past airfoils, as we have seen, 
modeled equations of turbulent flow have given evidence of their ability 
to capture the Hopf-type bifurcations typical of the onset of vortex- 
shedding. In application to the Inescapably three-dimensional flows 
typical of modern slender aircraft, modeled equations of turbulent flow 
will be called upon additionally to capture bifurcation phenomena such 
as the asymmetric vortex flows and vortex breakdowns that will appear 
as the parameters are varied over their extensive ranges. 
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FICURE CAPTION 


Fig. I Bonne- try conditions for loading due to • in angle of attack. 





